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Abstract. In recent high resolution N-body CDM simulations, it has been had found that nonsingular three- 
parameter models, e.g. the Einasto profile has a better performance better than the singular two-parameter models, 
e.g. the Navarro, Frenk and White in the fitting of a wide range of dark matter halos. A problem with this profile is 
that the surface mass density is non-analytical for general values of the Einasto index. Therefore, its other lensing 
properties have the same problem. We obtain an exact analytical expression for the surface mass density of the 
Einasto profile in terms of the Fox H-function for all values of the Einasto index. With the idea of facilitate the 
use of the Einasto profile in lensing studies, we calculate the surface mass density, deflection angle, lens equation, 
deflection potential, magnification, shear and critical curves of the Einasto profile in terms of the Meijer G-function 
for all rational values of the Einasto index. The Meijer G-function have been implemented in several commercial 
and open-source computer algebra systems, thus the use of the lensing properties of the Einasto profile in strong 
and weak lensing studies is straighforward. We also compare the Sersic and Einasto surface mass densities profiles 
and found differences between them. This impUes that the lensing properties are not equal for both profiles. 

1. Introduction 

The Cold Dark Matter (CDM) theory has become the standard theory of cosmological structural 
formation. Its variant the KCMD with ( Q^, ^a) = (0 -3, 0.7) seems to be in agreement with 
the observations on cluster-sized scales (IPrimackl 120031) . On galaxy/sub-galaxy scales has sev- 
eral problems, such as the discrepancy between observations and the results of numerical sim- 
ulations. The high resolution observations of rotation curves of low surface brightness (LSM) 
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20091). In contrast N-body 



CDM simulations predict a two para meter functional form for the density profiles with too high 
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densities (cusps) in the galatic center (iNavarro et al][l996 , 
ancy is called the cusp-core problem. 

Gravitational lensing is one of the most powerful tools in observational cosmolog y for probing the 
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regimes. In order to obtain the fit to the 



observational data in strong and weak lensing studies the most accurate density profile must be 
used. The first step before using a profile in lensing studies is to investigate the lensing properties 
of the profile. 



Recently N-body CDM simulations ( Navarro et al. 



Havashi & White 2008 



Stadel et al 



2009; 



Merritt et al 



2006 : 



Gao et al 



2008 



2004 

Navarro et al.l2010l) have found that the three-parameter 



profiles fit better to a wide rang e of dark mat ter halos. One of this profiles is the 
profile, a 3D version of the 2D 



Sersic 



Einastol(ll965h 



(119681) model used to described the surface brightness of 



galaxies. The Sersic profile can be written as: 



(1) 



where R is the distance in the sky plane, n the Sersic index, Y the mass-to-light ratio, 1^ the luminos- 
ity density at the effective radius R,,, b„ is a function of n that can be determine d from the condition 
that th e luminosity inside Re equals half of the total luminosity, for example IPrugniel & Simien 
J 19971) found b„^2n- 0.3333 + 0.009876/n. 
The Einasto profile is given by; 



p (r) = pe exp -da 



1 



(2) 



where r is the spatial radius, a is the Einasto index that determines the shape of the profile, da is 
a function of a which allows the calculation of the density pE inside an effective radius r^. In the 
context of dark matter halos this can expressed as: 



P(r) 



-■ p-2 exp — 
a 



— I -1 

r-2 



(3) 



where p_2 and r_2 are the density and radius at which p(r) cc r ^. Both radius and d ensities are 
related by p_2 - Pe exp(2/a - da) and r_2 - {adalTf. First iNavarro et al.l ( 120041) found that 



for haloes with masses fr om dwarfs to c l usters 0.12 < c < 0.22 with an average value of o' = 0.17. 

(l2008h found that there are trends for ue to increase with 



Havashi & White! (I2008h and 



Gao et al. 



mass and redshift, with g ~ 0.17 for ga l axy and g ~ 0.23 for cluster-sized haloes in the Millennium 



Simulation (MS) dSpringel et al.ll2005h. 



Navarro et al 



( 20101) f ound similar resu lts for galaxy-sized 
horn found that a ~ 0.3 



Gao et al 



haloes in the Aquarius simulation (iSpringel et al Also. 
for the most massive haloes of the MS. 

A problem with the Einasto profile is that in order to study its lensing properties numerical methods 
must be used because analytical expressions for this properties are not available. A semi-analytical 



Dhar & WilUams 



approx imation for the projecte d 2D projection of the Einasto was obtained by 
(I2OIOI) . A recent work done bv lBaes & Gentild (1201 lb demonstrated that is possible to write ana- 
lytical expressions for the deprojected Sersic model in terms of the Meijer G-function using MeUin 
integral transforms. 

In this paper, we present analytical expressions for the lensing properties of the Einasto profile. We 
apply the Mellin-transform method to derive an analytical expression for the projected surface mass 
density, the deflection angle, the lens equation, the deflection potential, the shear, the tangential 
and critical curves for all values of the Einasto index a in terms of the Meijer G-function. This 
function can be automatically handed by numerical routines implemented in computer algebra 
systems (CAS) such as the commercial Mathematica® and Maple® and the free open-source 
Sage and in the Python library mpmath. 
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This paper is organized as follows. In Section|2]we derive the surface mass density of the Einasto 
profile in terms of the Fox H-function and the Meijer G-function for values « = ^ and a - 
with n integer, and rational values of the Einasto index. In Section |3] we evaluate the total mass 
enclosed by this class of models using the density profile and the surface mass density obtained in 
the previous section. In Section|4]we use the result for the projected surface mass density to calcu- 
late the deflection angle, the lens equation and deflection potential for a spherically symmetric lens 
described by the Einasto profile in terms of the Meijer G-function. In Section|5]we derive expres- 
sions for the magnification, shear and the critical curves of the Einasto profile. We summary our 
conclusions in Section|6] We give an brief description of the Mellin transform-method in Appendix 
lAl In Appendix IB] we formulate all the properties of the Fox H-functions and Meijer G-functions 
that are used in this work. 



2. Analytical expression for the surface mass density of the Einasto profile 



The projected surface mass density of a spherically symmetric lens is given by integrating along 
the line of sight the 3D density profile: 



\J —CO 



p (^, r) dz. 



(4) 



where ^ is the radius measure fro m the centre of the lens and r - -\/^^Tz^. This expression can 
also be written as an Abel integral ( Binnev & Tremain3ll987 ): 

p(r) rdr 



By inserting equation Q into the above expression 
exp(^)st/s 



(5) 



S (x) - 2p_2?"-2e 



(6) 



having introduced the quantities x - ^/r_2 and s - r/r_2- 

The integral (|6]l can not be expressed in ter ms of ordinary functions for all the values of ag . 



However, using the Mellin-transform method (IMarichev 



1982 



Adamchick 



1996; 



Fikioris 



2007) 



is possible the exact calculation of one dimensional definite integrals. The most powerful feature 
of this method is that the result is a Mellin-Barnes integral. This integral for a certain combination 
of coeflicients is the integral representation of a Fox H-function or a Meijer G-function (for details 
see Appendix lAli . 

Using the Mellin-transform method with the integral (|6l), with z- I and the functions: 



f{s) = 2p_2'--2e" exp 
1 



-2s" 



Q < s < X 



-1 



elsewhere 



(7) 



(8) 



and its Mellin transforms: 



{M/)(M)-2p.2r_2e"a-'|- 



-u/a 



(9) 



Combining equations ( fTOt and (IA.4I ) with m = 2y and m = 1 /a yields: 



(10) 



Z (x) = V^p_2?--2e^'" a; 



2nij^ rH+y) 



(11) 



Comparing the last equation with ( IB.l II ) is possible to obtain an analytical expression in terms of 
the Fox H-function for the surface mass density of the Einasto profile: 



2(x) = ^p^^^-^s" xH 



(12) 



(1, 1) 

Writting the surface mass density as a Fox H-function has an inconvenient. The Fox H-function 
despite having a great potential for analytical work in Mathematics, sciences and engineering no 
numerical routines has been implemented yet. We prefer to describe the lensing properties of the 
Einasto profile in terms of analytical functions that have numerical routines akeady implemented 
to facihtate its use in strong and weak lensing studies. 



2.1. Surface mass density of the Einasto profile as a Meijer G-f unction 

The Meijer G-function meets the requirement pointed out before. A list of the relevant properties of 
the Meijer G-function can be found in Appendix|B] We can use this function to write expressions in 
analytical form for most of the lensing properties of the Einasto profile. The Meijer G-function had 
been implemented in several commercial and free available CAS. This means that using the Meijer 
G-function in lensing studies is just as simple as use other special functions Hke Hypergeometric, 
Gamma and Bessel functions for example. 



Using a similar procedure to the one used bv lBaes & Gentilel(l201 ll) to obtain an analytical expres 



sion for the luminosity density in terms of the Meijer G-function for all rational values of the Sersic 
index we proceed to do the same to derive an expression for the surface mass density of the Einasto 
profile for all values of the Einasto index. 



2.1 .1 . Einasto index witli values o' = ^ and c = f witli n integer 

The equation (fTTT i can be written in terms of the Meijer G-function for the Einasto index with 
values a = ^ and a - ^ with n integer. But first, one substitution is required using the Gauss 
Multiplication formula ( Abramowitz & Stegunlll970 ): 



A'-l 



nr(z' + ^) = (2;r)¥A^|-^^'r(A^z'), 

./=o 

with z' = z/N, N = 2m and z = 2my, we get: 

2m- 1 

Fd +2m3;) = (2m)5+2m, (27r)^'«r(l +3;) F(^ +3;] 
Substituting the last equation into (fTTT i. we obtain: 

2m— 1 



(13) 



(14) 



(2?:)'"-' 



'2bI^(-r-')nr(^^.)M-''.v 



(15) 



Comparing with the integral representation of the Meijer G-function (IB. lb we found an analytical 
expression for the surface mass density of the Einasto profile: 



where b is a vector of size | given by: 



(16) 



(17) 



This result indicates that the form of the surface mass density of the Einasto profile differs from the 
surface mass density of the Sersic model (equationlTJ in functional form. 



2.1.2. Einasto index with rational values 

Also is possible to write this expressions for all rational values of the Einasto index. Using m = p/q 
with p and q both integer numbers equation (fTTT i becomes: 



1 r^r(-i + ^3')r(i+2;,3') 

2 (x) = p-2r-2 ine ^ J — 



dy 



(18) 



T{l+qy) 

Substituting the three Gamma functions in equation ( fTSl l using the equation ( fT4l l. we obtain an 
integral and compare it with the definition of the Meijer G-function, we find: 



P-2r-2e 1 



(2n} 



P f^2p+q-\.Q 
q q-l,2p+q-l 



a 




b 


q^P 



where a and b are vectors of size q - I and 2p + q - 1 respectively given by: 

12 q-1' 
q'q''"' q 

2q-3' 



a 



^ {2/7' 2p 



2p - 1 113 
2q' 2q' 2q' 



(19) 



(20) 



(21) 



2 2q 2q 2q 2q 

It is immediate to verify that the equation (1% is equivalent to the equation ( fTSl l for Einasto index 
with values oe - \ - Using the properties of the Meijer G-function ( IB.3IIB.4b and ( IB. 61 ) is possible 
to demonstrate that equations (fT9T l and ( fTSl l are equal for Einasto index with values oe - \- 



2.1 .3. Simple cases: qe = 1 and oe -2 

In order to check our results, we projected the Einasto profile in the cases o- = 1 and a - 2 and 
compared the results with equation (fTSt . In both cases the values of a are outside the range favored 
by the N-body CDM simulations, but are practical to check the consistency of our calculations. 
For the case a = 1 we have: 



P(r) =p-2exp -2 — - 1 
\ [r-2 

Calculating the projected surface mass density using the equation (|5]l, we find: 

2 (x) - 2p_2''-2e^ X Ki (2x) 



(22) 



(23) 



where Ki (x) is the modified Bessel of the second kind. 
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Setting a = 1 in equation (fTSI l. we have: 



S (x) = p-ir-ie^ X Gq^j 



i _i 

2' 2 



(24) 





















v 







Substituting the equation (IB. 10b into (l24b we obtain the equation ( |23]) . 
In a similar way with the case a -2: 

p (r) = p_2 exp 

The projected surface mass density can be found using the equation (|5]): 

I(x)= V^p_2r-2e--"^' 

Setting a = 2 in equation (fTST l. we have: 



(25) 



(26) 



Z(x) = -v9rp_2'"-2e' -^Gq'j 



(27) 



Using (IB.9I I in the last equation this one reduces to equation i 
It is interesting to compare these two cases with the surface mass density of the Sersic profile 
2s (R) with the same values for the Sersic index 1 jm that for the Einasto index a. We also include 
the cases a - 0.5 {m - 2) and a - 0.2 (5) for the comparison. The Figure [T] shows 2s (R) for 
four values of m and Figure |2] displays £(x) for four values of a. In both it can be seen clearly 
that the respective index is very important in determining the overall behavior of the curves. The 
Sersic profile is characterized by a more steeper central core and extended external wing for larger 
values of the Sersic index m. For low values of m the central core is more flat and the external 
wing is sharply truncated. The Einasto profile has a similar behavior, with the difference that the 
external wings are most spread out. Also in the inner region for both profiles with low values of 
the respectively index we obtain larger values of and S. However, the Einasto profile seems to 
be less sensitive to the value of the surface mass density for a given a and radius and in the inner 
region than the Sersic profile. It is in this region where the lensing effect is more important and the 
difference in the surface mass density determines the lensing properties of the respectively profiles. 
Given this difference, we see that the lensing properties of the Sersic and Ei nasto profile are not 
equal. Studies of the lensin g properties of the Sersic profile had been done by 
Ehasdottir & MolleJ J2007[) . 



Cardonel(l2004 and 



3. The total mass enclosed 

The total mass enclosed in a halo described by the Einasto profile can be found by; 

M,„, =4;r I p(r)r^dr 
Jo 

Combining equations (O and (l28l l. we get: 
This result was also obtained by 



(28) 



(29) 



Cardone et al 



(l2005h . 



The Sersic Profile 
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Fig. 1. Sersic profile where Y/^ and Re are held fixed for four values of the Sersic index m. 
We can get the same result calculating the total mass projected on the sky plane: 

Jo 



Inserting equation (fT9] l into ( l30b , we find: 



a 


x' 


b 


q^P 



— -1 / 



(30) 



(31) 



Integrating the last equation using the formula ( IB. 71 ) for indefinite integration of the Meijer G- 
function: 



nor" e'i , 3 n'T"' r(f + Z;,) 



(32) 



We can write both products appearing in the numerator and denominator in equation ( l32l l using the 
Gauss multiplication formula ( fT3T l respectively as: 



n^fo'r(4?)n-^r(iii) 



2q 



2-1 V ? ; 



(33) 
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The Projected Einasto Profile 



10^ 



100 



0.01 



~ ~ — ^Ns*^^ 

^^^^ 


1 , , , 1 , , , 1 , 


1 , , , 1 

- 




^ " ■ " - . 






\ 

\ 


0.2 


- 








\ ^ 












\ 


0.5 




\ 






\ 






\ 


1 




V 




\ 






\ 






s 






\ 






\ 






\ 






a=2 \ 


\ 






\ 

\ 

\ 






\ 






\ 






\ 






\ 



20 40 60 80 100 

Radius 

Fig. 2. Projected Einasto profile where p-2r-2 and r_2 are held fixed for four values of the Einasto 
index a. 



Substituting equations (l33T l and (l34l i into (|32] | with a = ^ we obtain the same result $29[ for the 
total mass enclosed of the Einasto profile. This confirms that our calculations for the surface mass 
density of the Einasto profile are correct. 



4. The deflection angle, lens equation and the lensing potential 

4.1. The deflection angle and lens equation: 

In the thin lens approximation, the lens equation for a axially symmetric lens is: 
t^^^^-Dls& (35) 

where the quantities rj and ^ are the physical positions of the of a source in the source plane and 
an image in the image plane, respectively, a is the deflection angle, and D^, Ds and Dls are the 
angular distances from observer to lens, observer to source, and lens to source, respectively. 
With the dimensionless positions y = Dirj/Ds and x - ^/r_2, and dimensionless a - 
DiDisajDs^ the lens equation reduces to: 



y - X - a{x) 



(36) 



The deflection angle for a spherical symmetric lens is (iSchneider et alJll992l) : 



a(x) — — X cix — xK(x)ax 

X Jo ^cr/f X Jo 



(37) 



where k{x) = S (x) /Xcw/ is the convergence and Y^cru is the critical surface mass density defined by: 



c^Ds 



where c is the speed of light, G is the gravitational constant. 

Inserting equation (fT9] l into ( |37] |. we find the deflection angle of the Einasto profile: 



(38) 



a(x) 



P-2r-2e " 



- x^ (J^P+^-l.! 





^2./ 


b,-f 





(39) 



Introducing the central convergence, Kc, a parameter that determinate the lensing properties of the 
Einasto profile, defined by: 



Kc = 



^(x^Q) P-ir-ie 



p-2r-2e"r(ijpji- 



(40) 



and use it to write a (x) in terms of Kc'. 



a {x) — 



{2n)"-'T[^q^q\q 



£-1 



2 ^2;;+?- 1, 1 
q,2p+q 



b,-f 



2p 



(41) 



For Einasto index with values a - \ and a —\ with « integer, the last equation can be written as: 

1-1 



a (x) 



^2^il 



(27r)Hr(i) vsl 
The lens equation for the Einasto profile is then: 





x^ 







(42) 



(2n)P-'T{!^q^q\q 
which can be simplified to: 



q,2p+q 





x^i 


L 2q 


q^P 



(43) 



(27r)Hr(i) V^l 









x2 











(44) 



for Einasto index with values a - ^ and a - ^ with n integer. 

For a spherically sym metric lens being capab le of forming multiple images of the source a sufficient 
condition is Kc > \ dSchneider et al.lll992l) . In the case Kc < \ only one image of the source is 
formed. In addition to the condition Kc > I multiples images are produced only if | y |< ycrii 



dLi & Ostriker 



20021) . where jcrit is the the maximum value of y when x < or the minimum for 



X > 0. For singular profiles such as the NEW profile, the central convergence always is divergent, 
hence the condition Kc > 1 is always met, this implies that the NEW profile is capable of forming 
multiple images for any mass. Nonsingular profiles such as the Einasto profile are not capable of 
forming multiple images for any mass. Instead, the condition Kc > 1 sets a threshold for the lens 
mass required to form multiple images. 
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4.2. The deflection potential 

The deflection potential ^ (x) for spherically symmetric lens is given by; 



a {x) — 



dilr 
dx 



(45) 



We see from equation i45[ that can find the lensing potential simply integrating the deflection 
angle: 



t//(x) - I a{x')dx' 
Jo 

Inserting the equation (HTI) into ( |46] | and using the identity iB.ll . we found: 



2(27rf-'r(£)^2 y q\q 
which can be reduced to : 



1 - ^ 1 - ^ a 

2,' 2? 



„2? 



ip 







2' 2'^ 


x2 




2,^+2 


L 2' 2 





2(27r)"-'r(i) 

for Einasto index with values a - - and a = - with n integer. 



(46) 



(47) 



(48) 



5. Magnification, shear and tlie critical curves 

The gravitational lensing effect preservers the surface brightness but causes variati ons in the shape 



and th e solid angle of the source. Thereby, the source luminosity is amplified by (ISchneider et al 



1992) 



1 



(49) 



where k {x) is the convergence and y (x) is the shear The amplification has two contributions one 
from the convergence which describes an isotropic focusing of light rays in the lens plane and the 
other is an anisotropic focusing caused by the tidal gravitational forces acting o n the light rays, 
described by the shear. For a spherical symmetric lens, the shear is given by (IMiralda-Escude 



199 ih : 



t{x)-W _ 

7 (X) = ^ K- K 



where 



2 r"" 

Z(x) = — x'I.(x')dx' 
Jo 



(50) 



(51) 



is the average surface mass density within x. 

The magnification of the Einasto profile can be found combining equations ( fT9] l, ( |49] l, (ISOl l and 
dlB- We get: 



where 



k(x) = 



a 


^2, 


b 





(52) 



(53) 
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The last equations reduce to: 







L ' 2q 





(54) 



k{x) 



k(x) = 



(27r)Hr(i) V^U 



X G ° 



jc G° 



a 








b 







x2 







(55) 



(56) 



for Einasto index with values a - ^ and a - with n integer. 

The magnification may be divergent for some image positions. The loci of the divergent magnifica- 
tion in the image plane are called the critical curves. For the Einasto profile we see from equation 
(|52] | that has one pair of critical curves. The first curve 1 - = is the tangential critical curve 
which correspond to an Einstein Ring with a radius called the Einstein radius. The second curve 
I +'k-2K - Qis the radial critical curve which also defines a ring and its correspond radius. In both 
cases the equations must be solved numerically. 



6. Summary and Conclusions 

In this paper, we have derived an analytical expression for the surface mass density of the Einasto 
profile using the Mellin transformed. This expression can be written in terms of the Fox H-function 
for general values of the Einasto index a. The same expression can be written in terms of the Meijer 
G-function for all rational values of the Einasto index, with a simphfication for values a = ^ and 
a - with n integer of the Einasto index. One we obtained an analytical expression for the surface 
mass density we also derived in terms of the Meijer G-function other lensing properties: deflection 
angle, lens equation, deflection potential, magnification, shear and critical curves of the Einasto 
profile for all rational values of the Einasto index, with a simplification for values a - ^ and a - ^ 
with n integer of the Einasto index. Our analytical results can be used to investigate further the 
lensing properties of the Einasto profile taking advantage of the fact that the Meijer G-function is 
a very well studied function in the literature. 

We compared the Sersic and Einasto surface mass density profiles using the equivalent values for 
the Sersic m and Einasto a indexes and where the quantities T/^, Rg and p-ir-i, r^2 are held fixed. 
We found that both profiles have similar behavior determined by the index value. However, we 
noted that for the Einasto profile the external wings are most spread out and seems to be less 
sensitive to the value of the surface mass density for a given Einasto index and radius in the inner 
region than the Sersic profile. This feature is key because it is in this region where the lensing 
effect is more important and the difference of the surface mass densities implies a difference in the 
lensing properties of the two profiles. 

Our results can be used in strong and weak lensing studies of galaxies and clusters where dark 
matter is to believed the main mass component and the mass distribution can be assumed to be given 
by the Einasto profile. The implementation of this results is easy because the Meijer G-function is 
available in several commercial and open-source CAS. The performance of this nonsingular three- 
parameter model in fitting the 3D spatial densities in high resolution N-body CDM simulations is 
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better than the singular two-parameter NFW profile makes very promising its use in strong and 
weak lensing studies. The constant increasing computational power available opens the possibihty 
of using most realistic and sophisticated profiles hke the Einasto profile for lensing studies and 
marks a route to obtain a satisfactory solution to the cusp-core problem. 
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Appendix A: The Meijer transform-method 



The Mellin transform-method dMarichev 



1982 : 



Adamchick 



1996 : 



Fikioris 



2007h uses the MeUin 



integral transform for the integral evaluation. 

The Mellin transform of a function / (z) is an integral transform defined by: 



{Mf]{ 



«)= z"- 



'f(z)dz 



(A.l) 



if the integral exits. 

It is clear from the definition that the Mellin transform does not exist for all functions such as the 
polynomials, the integral does not converge. The Mellin transform when it does exits it converges 
in a vertical strip in the complex z-plane. This strip is called the strip of analyticity (SOA). 
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The inverse Mellin transform is defined by: 
f(^^^^ij z-"{Mf}(u)du 



(A.2) 



where the contour of integration C is a vertical line in the complex z-plane and must be placed in 
the SOAof/(z). 

Given two functions / (z) and g (z) the Mellin convolution is defined by: 



(A.3) 



It is well know that the Laplace or Fourier transform of the product of two different functions is the 
convolution of the respectively transform. In the case of the Mellin transform we have: 



f(y) 



(A.4) 



if z = 1 this formula is know as the Parseval's theorem for the Mellin transform. 
The most important feature of the Mellin transform-method is that using the equation (IA.4b inte- 
grals of the type, 



(A.5) 



can be written as an inverse Mellin transform. With the requirement that / and g should be of the 
hypergeometric type and consequently their Melhn transforms can be written as products with the 
form r (fl -H Au) or [F (a + Am)]"' with the A's being real numbers, the resulting integrals are of the 
Mellin-Barnes type and then can be written in terms of the Fox H-function for A 1 or the Meijer 
G-function for A = 1 (see AppendixlBl). 



Appendix B: The Meijer G function and its properties 



The Meijer G-function is a very general, analytical function introduced by iMeiien (Il936h which 
includes most of the spec ial functions as specific cases. It is defined in terms of the inverse Mellin 
transform ( ErdelvJl954 ) by: 



a 






fll. 


..,ap 




z 




z 


b 






■ ■,b, 





1 r n';Lir(^ + .)n;iir(i 

r[aj + s) 



^ A J=m+l \ J / A 1 j=n+\ 



z "ds 



(B.l) 



where C is a contour in the complex plane, T{s) is the Gamma function and a and b are vectors of 
dimension p and q, respectively. 

The basic properties of the Meijer G-function are too numerous to be mention here. We only pro- 
vide a short list of the most relevant properties for this work. 

A Meijer G-function with p > can be transformed to another G-function with p < q: 



\-bu-A-bq 

1 — fll, I — Op 



(B.l) 
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Other property is that if one the parameters of a and b appears in both the numerator and denomi- 
nator of the integrand, the order of the Meijer G-function may decrease and the fraction simpHfied. 
The positions of the parameters dictates which order m or n will decrease. For example if - bj 
for some k - 1,2, n and j = m + l,m + 2, q , the orders p, q and n of the Meijer G-function 
will decrease: 



01,02, ...,ap 
bi,...bq-],ai 



1 

p-\,q-\ 



Cl2, Clp 
bl,...,bq-l 



(B.3) 



In the other case if at - bj for some k — n + l,n + 2, ...,p and j - 1,2, ...,m, the orders p,q and m 
of the Meijer G-function will decrease: 



b\,b2, -ba 



= G 



■m-l,n 
p-\,q-\ 



a\, ...,ap-i 
b2, ■■■,ba 



The order reduction formula for the Meijer G-function is: 



(B.4) 
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■ ■,b. 
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{2nf-'^' 



■^km,kn 
kp,kq 



a\/k, (ai + k - 1) //:, ap/k, (^Gp -\- k - /k 

bjk,...,(bl+k-l)/k,...,b,/k,...,{bg+k-l)/k 



The multiplication by powers of z is another property: 



p,q 



Cl\, Gp 

b\,...,bq 



a\ + a, Gp + a 
b\ + a, bq + a 



(B.5) 



(B.6) 



Among the indefinite and definite integrals of the Meijer G-function one has the following: 

1 — a,a\, a 



^P,i 



ai. 
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az 


bi. 


■ ■,b. 





z dz-z G„„ 



(B.7) 



Jo '-'P'1 



bi. 



z"-'dz 



\{l,T{bj + a)Wj=j[l-aj-a) 



(B.8) 



n^.„,,,r(i-/,,-a)n;.„„r(a, + a) 

A short list of relations between the Meijer G-function and some elementary and special functions 



-,1,0 
^0,1 



-.2,0 



exp (-z) z 



bi,b2 

A more complete list can found in 



(B.9) 



(B.IO) 



Bateman & Erdelvi 



(119531) and the Wolfram Functions Sitel^ 



http://functions.wolfram.com/HypergeometricFunctions/MeijerG/ 
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The Fox H-function is a generalization of the Meijer G-function introduced by 
defined in terms of an Mellin inverse transform: 



Foxl (1196 ll) . It is 



Tjm.n 



(a, A) 
(b,B) 



_ Tjm.n 



(fli, Ai), ...,(flp, Ap) 
(bu Bi),...,ib,, B,) 



J_ r Y\%T{bj + Bjs)Y\y_,Y{\-aj-Ajs) 
X n^.,„,, r(l - bj - Bjs) U%„,, r{aj+Ajsy 



(B.ll) 



Srivastavaet al.ll 1982 



Hai & Yakubovich 1991 



Kilbas & Saigoll2004t 



More properties and apphcations of the Fox H-function can be found in ( Mathai & Saxena 

Mathai et al.ll2009l) . 
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